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The contributions of different color flows into the gluonic Sc-nieson production 
cross section has been calculated. This study is essential to simulate Bc-raeson pro- 
duction with the help of Pythia program. The essence of matter is that in the frame 
■ work of the Lund model used by Pythia the hadronization way of the final par- 
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O ■ tons and hadronic remnants depends on color flow type. The modified method for 
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calculation of the color flow contributions has been proposed. 



1. INTRODUCTION 



A simulation of particle production processes at modern hadronic experiments is practically 
impossible without Monte-Calro methods. The fact of the matter is that the integro-differential 



cn 
o 

^ . equations, which are used to describe the evolution of initial partons and the hadronization of 

^ parton interaction products, are very complicated. Recent time PYTHIA [1] is the most dependable 

P_^' program for Monte-Carlo simulation of a particle production at high energies. This software 

! simulate all three studies, on which the hadronic production can be conventionally subdivided: the 



initial parton evolution, the hard subprocess of the initial partonic interaction, the hadronization 
of the hard interaction products and hadronic remnants. The hard subprocess is calculated in the 
frame work of perturbative QCD. The hadronization in Phythia is described by Lund model [2]. 
In the frame work of this model one supposes that outgoing color charges stretch the connecting 
force fields (color strings), and it is these fields that eventually break up to produce the final 
state hadrons. The force field is always stretched from a color triplet to a color antitriplet. Color 
octets (for example, gluons) are treated as an excitation of the string. A technique which allows 
to calculate the particular color flow contribution into the cross section has been performed in [2]. 
Also it has been demonstrated in [2] that at infinite color number limit Nc —>■ oo the interference 
terms between different color flows equal to zero. 

A wide set of standard hard processes can be simulated by PYTHIA, however this set does 
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not contain all processes, which would be interesting to study at high energies. That is why in 
the recent versions of PYTHIA a user have a possibility to include his own matrix element into 
the PYHTIA program. Therefore the problem of color flow separation appears. It is important 
to know the contribution of different color flows into the cross section, because each color flow 
corresponds to its sole hadronization way. 

In this work the color flows for the process of gluonic production of Sc-meson have been studied. 
The modified calculation method of the color flow contributions has been proposed, which is slightly 
differ from the traditional one. 

It is worth to note that it is somewhat difficult to calculate the cross section for the process 
of gluonic i^c-meson production, because Bc-meson production mechanism can not be treated as 
fo5-pair production followed by 6-quark hadronization into Bc-meson. Therefore Bc-meson pro- 
duction can not be described by the fragmentation function b ^ Be (the fragmentation function 
calculations have been performed in [3]). It has been shown that under conditions of present-day 
and planed experiments the recombination mechanism of Sg-mesons production dominates. To 
estimate the recombination contribution into the cross section one need to calculate 36 tree-level 
Feynman diagrams of order 0(0;^). These calculations have been done independently by several 
research groups [4-9]. The results of [4-7] are in a good agreement with each other. 

An interest in Sp-meson production in hadronic interactions is increasing now due to the fact 
that LHC will have started to work in the nearest future. The last versions of PHYTIA based 
generator SIMUB [10] include the useful subroutines to simulate Be production, which are based 
on codes of research groups [4, 5]. The fruitful scientific contacts between the author and SIMUB 
developers have caused the presented study. 

2. THE SEPARATION OF COLOR FLOWS 

The Feynman diagrams for the process gg ^ Bc-\- X are shown in Fig. 1. It is useful to treat 
the diagram with four gluonic vertex as three different diagrams, because that diagram contains 
three different color structure. The color parts of the diagrams are given by following equations: 

rp fniS2Tl2 fn3n2giJ.niJ.n3!: _ 

J-1 — J J ^CC^feb "&C) 

rp fnig2n2 fn3n2gij.n3j.n1r 

rp _ i-nigig2 fn2n1n3.Ln2j.ni _ 
^3 — J J i^a. I'cc "be, 



rji^ _ yniS2n2 yn2n35i^"_3^ni^^_ 
rp^ _ jg2nin2 jn2nsgijJii^nz 
rp^ _ jnin2n3 yn3925i^n_i^n2^^_ 

-'IS — ''bli''hVcl2''l2c"bc, 

J 16 - hlJ'l{b^cl2^l2<Pbc, 

T — +91 +ni J.92 j.ni r 
^17 - hh^iib^cl2h2C^bc, 

_ 7-91 T-ni 7-92 j.ni r 
-'is - hh^l{b^cl2^l2C°bc, 

rp j.ni J.52 j.gi j.ni r 

J 19 - hh'^lib''cl2h2C^bc, 

_ j.ni ,§2 j.ni ,gi sr 
J 20 - hh^lib^ch'^hc^bc, 

J 21 - hh'^i{b^cl2h2<Pbc, 

T — +92 j.ni j.ni ,gi r 
J 22 - Viib^d2h2<Pbc, 

J 23 - hhhih^l^b^cc^bc, 

rp _ J.91 j.ni 1 92 T-ftiX 
J 24 - hhhih^l^V'ccObc, 

T _ +92 ^ni ^ni r 
J 25 - hhhih^i^b'^ccObc, 

T _ +'^1 +92 ^31 ^ni c 
J 26 - hh\l2^l^b^cc(^bc, 

T _ +52 ^ni .91 .ni r 
-^27 — hhhih^l^l^cc'Jbc, 

T _ +92 ^91 +"1 +rai X 
J 28 - hhhih'^l^b'^cc'^bc, 
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^29 


_ fni.gi .92 fmc 

— ''Ch''hl2''l2c''bb"^''' 


^30 


_ j.gi A.ni ,g2 ^ni r 
— ^di^hh^hc^bb^^^-' 


7^31 


_ 4-91 +92 r 
— ^ch^hh^hc^bb"^' 


T32 


— ^dlhll2h2C^bb"^^' 


^33 


_ f92 -r»i -31 j.ni r 
— ^chhibhic^b 


^34 


_ J.g2 j.gi j.ni ,ni r 
— ^dlhll2h2c''bb"''^' 


^35 = 


■ fnigig2fn2 +"if"2X 


336 — 


^ fn1g1g2j.n1 ,n2 J.n2j: 

'^bh^l-^b'^cc^bc' 


T37 — 


Afni9ig24.n2 j.ni ,n2 r 

^cli^hc^bb^bc' 


Tss — 


Afn1g1g24.n1 ,n2 +"2^ 



where upper indexes (71, (72 are color states of the initial gluons, lower indexes 6, 5, c, c are color 
states of 6-, 6-, c-, c-quarks correspondingly and 5bc is a color wave function of the 5c-ineson ( the 
normalization coefficient l/\/3 is not written for the sake of simplicity). 

Let us study, for example, the color part of diagram 1 (see Fig. 1. In this diagram the initial 
guons exchange a gluon in t-channel and split into the quark-antiquark pairs): 

_ jnig2n2 yn3n2gi^ni^n_3^^_ _ jnig2n2 Jra3n29i ^^"i^"3^_ 

Using the identity tH^ — tH"' — if^^H"^, one can found: 

jnig2n2 yn3"2Si^ni^r!.3 _ _^^52^n2 _ ^02 ^52 ^^^"2^51 _ ^fli^ni^ _ 



.1^52^31 _ ^1^2^91 _ i^</2^9l _ (1^9192 _ 1^92^91) 

3 6 6 M 6 ^ 



Thus: 



^ _ 1(59192 _ ^^92^91, 
4 2 
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Q 1 
TL — +924.91 _ 

The first term (see the scheme (2) in Fig. 2) corresponds to the case, where a color of the 
second gluon [g2 flows to c-quark, an anticolor of the first gluon [gi) flows to 6-quark, and a color 
of the first gluon annihilates with an anticolor of the second gluon (sum over k). The second term 
corresponds to the case, where colors and anticolors of the initial gluons annihilate, and a color and 
an anticolor of c-quark and 6-quark are produced from a vacimm (see the scheme (3) in Fig. 2) . In 
addition to the terms described above, the term tcfc^S contributes to another diagrams. The latter 
one corresponds to the case, where a color of the first gluon (^fi flows to c-quark, an anticolor of 
the second gluon {g2) flows to 6-quark, and an anticolor of the first gluon annihilates with an color 
of the second gluon (see the scheme (1) in Fig. 2). There are no color flows for discussed process 
but such as three ones described above. This color flow separation has been done by following the 
recipe given in papers [2]. 

Nevertheless, it would be better to base on more fundamental QCD principles to describe color 
flows. Let us consider the term which corresponds to the production of c-quark and h- 

quark in a color singlet. Naturally, a color sting stretch between these two quarks. It is worth to 
note that the term t^t^j^i contains the singlet part too: 

tH^ = l^"^ + i(d«''- + ip^y. (1) 

Therefore "one part" of the singlet is hadronized in one manner, and "other part" of the singlet 
is hadronized in other manner. We think that it would be better to treat the total color singlet 
contribution as the separate color flow. Two other color flows would be composed of two color 
octet states d and /. 

That is why we redefine the color flows as follows: 

1. A color of the first gluon flows to c-quark, an anticolor of the second gluon flows to 6-quark, 
an anticolor of the first gluon and a color of the second one annihilate: 

2. A color of the second gluon flows to c-quark, an anticolor of the first gluon flows to 6-quark, 
a color of the first gluon and an anticolor of the second one annihilate: 
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3. A color and an anticolor of the initial gluons annihilate, a color of c-quark and an anticolor 
of 6-quark are produced from vacuum: 

In our point of view these definitions of colors flows are more physically justified, because the 
singlet state contribution is separated from the octet contributions and there is no an interference 
term between the singlet color flow and other flows. 

The definitions described above slightly differ from the traditional ones: ^cfe^S' ^ck^ll ^^'^ ^^^^^^cl- 
However, an accuracy of the color flow separation is about 1/Nc, where Nc is a color number. For 
an arbitrary Nc the equation (1) looks like follows: 

tH'' ^-Ls''^+- {(T^'' + ir^^)^, (2) 
INf, 2 

and one can conclude that at 1/Nc oo, the both definition sets lead to the same results. 
Thus the color part of matrix element n can be given by the following expression: 

Tn = \{d^''''' + ir''^)t% ■ + - ir^^'')t% ■ S„ + . C„, (3) 

From (3) one can easy obtain the color matrix averaged over the initial color states and summed 
over the final ones: 

Mran = ^J^{D + F) ■ A^A^ + {D + F) ■ BmBn + {D - F) ■ {AmBr, + BmAn) + S ■ C„C„) , (4) 

where D — 5/3, F — 3, and S — 24. Vectors A„, B^ and C„ are performed at the table 1. 

In our approach the color matrix corresponded with an interference between the color flows 

{d<^^92k ^ ■p,g2k^^k_ ^^si32fc _ ^f9ig2k-^-i-k_ i^g^^g simple form: 

=^{{D-F). (A^Br, + B^A^)) (5) 

It is worth to note, that the formula (3) can be easily rewritten in the traditional color flow 
definitions: 

Tn = ■ A^ + t%t% . B^ + 5'^-^'^5ci ■ (C„ - ^l!L±^). (6) 

One can see that the expressions for the color flow (1) given by (3) and (6) differ from each 
other only in a common coefficient. It is easily to show that for the color flow (1) the ratio for the 
matrix element from (3) to one from (6) is 7 : 8, as well as for (2). 
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3. CALCULATION RESULTS 

The cross section distributions over exit angles of final particles {Be, b and c) have been shown 
in Fig. 3 and 4 for the different color flows and for the interference term at the gluon interaction 
energy 25 GeV. This energy value has been chosen because it is common value for ^c-meson 
production at LHC. It is worth to mention that intuitive ideas about color flows correspond to 
the calculation results. Indeed, it is clear from Fig. 3c and 4c that 6-quark moves mainly in the 
direction of gluon, which transfer an anticolor to 5-quark, Also one can see from Fig. 3f and 4f, 
that c-quark moves mainly in the direction of gluon, which transfer a color to c-quark. It is 
not unexpected, that the color flow (3), which corresponds to a singlet state, is symmetrically 
distributed over the angles. 

The interference term between flows (1) and (2) is small (see Fig. 3b, d, f H 4b, d, f). For 
the pseudoscalar i?c-nieson production at 25 GeV the interference contribution is negative for all 
values of the exit angles. For the vector meson production the interference contribution is negative 
too with the exception of peripheral regions of the angular distributions. The total interference 
contribution into the vector meson production is negative, as well as into the pseudoscalar meson 
production. Our calculations show that at low energies the interference contribution becomes pos- 
itive for the vector meson production and remains negative for the pseudoscalar meson production. 
An absolute value of the interference contribution is small at all interaction energies. 

As it was mentioned above, the contributions of the color flows (1) and (2) in our separation 
scheme differ from the traditionally determined contributions only by the common coefficient. It 
is not so for the color flow (3). The contribution of (3) in our approach one in traditional approach 
differs from each other in shape too. The distributions over exit angle of Sc-meson in the frame 
work of traditional and our approaches have been presented in Fig. 5. One can see that in our 
approach the contribution of (3) is larger than this contribution in the traditional scheme. The 
point is that in our approach there is no interference between the color flows (1) and (3), as well 
as between (2) and (3). Therefore the contributions, which correspond in traditional scheme to 
the interference between (1) and (3), as well as between (2) and (3), transfer to color flow (3). 

4. CONCLUSIONS 

The proposed calculation method minimize interference terms between color flows for the pro- 
cess gg Bc + b + c ^ Bc + X . That is why interference terms can be neglected in our approach. 
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Furthermore the method under discussion allows us to ascertain more clearly how color states of 
final partons correspond with color flows. The redistribution of interference terms leads to the 
amplification of the color flow (3) . The contribution of this color flow into the central kinematic 
region becomes comparable with the contributions of color flows (1) and (2). Thus haronization 
features of the outgoing partons b and c could be changed. 
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Table 1. The vectors An , Bn, Cn corresponded with the color flows (1), (2), (3) in Fig. 2. 



n 


An 


Bn 




1 





-3/2 


-1/2 


2 


-3/2 





-1/2 


3 


3/2 


-3/2 





4 





3/2 


1/2 


5 


-3/2 





-1/2 


6 


3/2 


-3/2 





7 








1/4 


8 





-3/2 


-1/4 


9 


3/2 





1/4 


10 








-1/4 


11 








-1/4 


12 





3/2 


1/4 


13 


-3/2 





-1/4 


14 








1/4 


15 





-1/6 


2/9 


16 





-1/6 


-1/36 


17 





-1/6 


-1/36 


18 





4/3 


2/9 


19 


4/3 





2/9 



n 




Bn 


Cn 


20 


-1/6 





-1/36 


21 


-1/6 





-1/36 


22 


-1/6 





2/9 


23 


4/3 





2/9 


24 


-1/6 





-1/36 


25 


-1/6 





2/9 


26 





4/3 


2/9 


27 





-1/6 


-1/36 


28 





-1/6 


2/9 


29 


-1/6 





2/9 


30 


-1/6 





-1/36 


31 


4/3 





2/9 


32 





-1/6 


2/9 


33 





-1/6 


-1/36 


34 





4/3 


2/9 


35 


4/3 


-4/3 





36 


-1/6 


1/6 





37 


-1/6 


1/6 





38 


4/3 


-4/3 
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Figure 2. Color flows for the process gg ^ Bc + c + b. Color strings stretch as follows for the color flow 
(1): [c-quark — > the remnant of hadron, which contained the gluon gi — > the remnant of hadron, which 
contained the gluon 52 — > 6-quark]; for the color flow (2): [c-quark — > the remnant of hadron, which 
contained the gluon 52 — > the remnant of hadron, which contained the gluon 51 — > 6-quark]; for the color 
flow (3): [c-quark — > 6-quark], [the remnant of hadron, which contained the gluon gi the remnant of 
hadron, which contained the gluon 52]- 



13 



da / d cos Q nb dcrint/c^cos 6^^, nb 




-0.5 0.5 1 -1 -0.5 0.5 1 

COS 6c COS 6c 



Figure 3. The cross section distributions for the pseudoscalar Bc-meson production at 25 GeV {gg 
Be + b -\- c) for the different color flows and the interference term: over the exit angle of Sc-meson ( a) and 
b) ) ; over the exit angle of 6-quark ( c) and d) ) ; over the exit angle of c-quark ( e) and f) ) . A solid curve in 
plots a), c) and e) denotes the color flow (1); a dashed curve denotes the color flow (2); a dotted one denotes 
the color flow (3). The interference contribution is shown in plots b), d) and f). 
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da / d cos Q nb dcrint/c^cos 6^^, nb 




-0.5 0.5 1 -1 -0.5 0.5 1 



COS 6c COS 6c 



Figure 4. The cross section distributions for the vector Bc-meson production at 25 GeV for the different 
color flows and the interference term. The designations are the same as in Fig. 3. 




Figure 5. The cross section distributions over the exit angle of pseudoscalar (a) and vector (b) Be- meson 
for the color flow (3): in our approach (soHd curve) and in the traditional one (dashed curve). 



